Quasiparticle transport equation with collision delay. I. Phenomenological approach 
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For a system of non-interacting electrons scattered by resonant levels of neutral impurities, we show 
that viriai and quasiparticle corrections have nearly equal magnitudes. We propose a modification of 
the Boltzmann equation that includes quasiparticle and viriai corrections and discuss their interplay 
on a dielectric function. 



I. INTRODUCTION 

Elastic scattering of electrons by impurities is the sim- 
plest but still very interesting dissipative mechanism in 
semiconductors. Its simplicity follows from the absence 
of the impurity dynamics, so that individual collisions are 
described by a motion of an electron in a fixed potential. 
On the other hand, due to a large variety of impurities 
and their accessible concentrations, impurity-controlled 
transport regimes span from simple response character- 
ized by a mean-free path to a weak localization. 

Let us recall a quasiclassical picture of impurity con- 
trolled transport. The basic effect of impurities on trans- 
port in crystals consists in abrupt changes of directions 
of electron trajectories. Within the Boltzmann equation 
(BE), this effect is described by scattering integrals. At 
higher concentrations, impurities influence a band struc- 
ture. This correction can be built into the BE if one ac- 
cepts that elementary excitations are not simple electrons 
but electron- like quasioarticles in the spirit of the Landau 
theory of Fermi liquids HI Finally, impurities attract/expel 
electrons to/from their vicinity what reduces/increases 
a density of freely traveling electrons. Such changes in 
effective density of electrons are covered by viriai cor- 
rections that are accounted for via non-local (in time 
and space) corrections to scattering integrals. Although 
quasiparticle and viriai corrections to the BE are known 
for decades, transport theory that would include both 
in the same time is still missing. Our aim is to fill this 
gap. In this paper we focus on an intuitive approach. 
In the second paper of this series we confirm equations 
presented here from the quantum statistics. 

To introduce basic concepts, we first discuss classi- 
cal viriai corrections to the BE, Sec. II. In particular, 
we show that non-locality of scattering events can be 
described in terms of a collision delay. In Sec. Ill we 
review Wigner's concept of collision delay and estimate 
magnitude of viriai corrections for resonant levels in III- 
V semiconductors. In Sec. IV we show that viriai cor- 
rections go hand-in-hand with quasiparticle corrections 
having nearly equal magnitudes. In Sec. V, an intuitive 
modification of the BE and of the most important observ- 



ables (density, current, energy) is proposed. In Sec. VI, 
we discuss viriai corrections to transport coefficients. In 
Sec. VII we summarize. In Appendix A, we show that 
non-self-consistent treatment used in this paper and more 
advanced self-consistent treatment are equivalent within 
assumed precision. In Appendix B, we derive the derived 
optical theorem which explains comparable magnitudes 
of the viriai and the quasiparticle corrections. In Ap- 
pendix C, we verify that the presented modification of 
the BE is consistent with the equation of continuity and 
the energy conservation. 

II. CLASSICAL COLLISION 

Since the quantum-mechanical theory with intuitively 
clear viriai corrections is still missing, the only experience 
for non-equilibrium systems one can gain from viriai cor- 
rections to the BE in the classical statistical theory of 
moderately dense gases. Within accuracy to the second 
order viriai coefficient, these corrections were introduced 
already on the break of centuries by Clausius.0 Here we 
modify his approach in two aspects. First, instead of bi- 
nary collisions of molecules we assume electron-impurity 
events. Second, instead of space non-locality, we refor- 
mulate viriai corrections in terms of time non-locality. 

According to Clausiusp one has to take into account 
that two colliding molecules arc not at the same space 
point, but at a distance of sum of their radii. In other 
words, the scattering integral has to be non-local in 
space. 

Similar argument about non-locality of scattering 
event applies to collisions of electrons with impurities. A 
sketch of a classical trajectory of a colliding electron is in 
Fig. [j]. Before the electron reaches the impurity potential 
of a finite range, its trajectory is a straight line. Then it 
makes a curve in the impurity potential and again follows 
a straight line in a new direction. Within the BE, this 
process is approximated by an effective event that is local 
in time and space. Of course, within the local approxima- 
tion, one has to sacrifice dynamics of the electron during 
the collision. More serious neglect follows from the fact 
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that within the local approximation the asymptotic mo- 
tion along outgoing line cannot be properly matched with 
the motion along the incoming line. 

Let us find a correct matching of incoming and outgo- 
ing lines. To this end, we extrapolate the incoming and 
outgoing lines and find their crossover X. In general, 
such a crossover need not exist, however, it always exists 
for spherical potentials to which we limit our attention. 
The crossover X gives us the coordinate at which we have 
to place the effective event. As one can see in Fig. [j], the 
crossover X does not coincide with the centre of impu- 
rity. Using the local scattering integral of the BE, the 
scattering event is placed in the center of impurity, thus 
the shift of the center of scattering event is the first ne- 
glect that influences motion of electron in the asymptotic 
region. 

The second neglect of the local approximation is not 
visible from the trajectory in Fig. |l[ The electron fol- 
lowing the extrapolated incoming line passes through the 
crossover X at time t m . The electron on the extrapolated 
outgoing line passes X at time t out . In general, t ln ^ t° ut , 
however, the time locality of scattering integrals in the 
BE means that t in = t out . 

One finds two misfits: in the position of the event, and 
in the matching of times. Although these two misfits usu- 
ally come together, their consequences can be discussed 
separately by using special models. 



A. Point traps 

First, we assume impurities of a negligible volume with 
a capability to trap electrons for a certain time. In this 
case, there is no misfit in the position but £ out — t m — 
A t > 0. 

Using intuitive arguments, the collision delay A t can 
be incorporated into the scattering integrals of the BE. 
A balance equation of the Boltzmann type for scattering 
by impurities reads 



df_ + k_dj_ _ d^_df_ 

dt m dr dr dk 



dp 
dp 



P pk f(p,r,t%) 
Pk P f (k,r,t™) , (1) 



where f(k,r,t) is a distribution function in the phase 
space, r is a coordinate, t is a time, k and p are mo- 
menta, Pfcp is the scattering rate from k to p. Since dis- 
tributions in the scattering integrals correspond to initial 
conditions, t™ p is a time at which an electron enters the 
scattering from k to p. 

In the scattering-out event [the second term on the 
r.h.s. of (]]])], an electron of momentum A; enters a col- 
lision at t k n leaving at with momentum p. The 



scattering-out integral gives a probability that at time 
t an electron leaves the momentum k. This happens at 
the beginning of the collision, thus t k n p 



In the scattering-in event [the first term on the r.h.s 
of (El)], an electron of momentum p enters a collision at 
t p \ leaving at with momentum k. The scattering- 
in integral gives a probability that at time t an electron 
enters the momentum k. This happens at the end of 
the collision, thus t°f = t. From t°f - t p \ = A t (p, k), 
one finds that t p \ = t — A t (p,k). The time argument 
in the scattering-in is thus shifted by the collision delay 
A t (p, k). A modified BE then reads 

df k df d<\>df [dp 

dp 



(2tt>- 



-P kp f(k,r,t). (2) 



Electrons trapped by impurities are excluded from free 
motion. With a finite collision delay, one has to deal with 
two distinguished local densities of electrons. Beside the 
physical density n = N/Cl (number of electron N per 
volume f2), there is an effective density 



n fTee {r, t) 



dk 



f(k,r,t) 



(3) 



which equals the local density in the free space between 
impurities. 

For finite collision delay A t , a share of electrons 
trapped by impurities can change in time. Accordingly, 
the free density nf rco does not conserve. From (||) one 
finds that in a homogeneous but non-stationary system 



dn {l , 
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dk dp 
(2tt) 3 (2tt) 3 
dk dp 



P pk (f(p,t-A t )-f(p,t)) 

OA KM 

(27r) 3 (27r) 3 pfe * dt 
d f dk dp 

The quantity that conserves is the physical density 
which differs from the free density by the density 



dk dp 
(2^)3 (27y 



■P pk A t f(p,t) 



(4) 



(•5) 



(6) 



that is correlated with impurity positions. 

Note that the scattering mechanism enters relation be- 
tween density n and distribution /. Without virial cor- 
rection (here represented by correlated density) , the func- 
tional n[f] is independent of scattering, since n — n.f ree . 
In the presence of virial corrections one has to keep in 
mind that a density of freely traveling electrons does not 
equal the physical density. 
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B. Hard spheres 

As the second example, we discuss hard-sphere impuri- 
ties. In this case, the incoming and outgoing lines have a 
crossover at the sphere surface. The times match exactly, 
^out _ £in rpj^ Qn jy misfit results from the fact that the 
crossover is not at the centre of the impurity but shifted 
by the sphere radius. Here we show that the crossover 
offset can be reformulated in terms of an effective time 
mismatch so that one can use unified description of col- 
lisions with point traps and with hard spheres. 

A collision with a hard sphere is schematically shown in 
Fig. ||. The real electron trajectory follows the full line. 
The scattering integral of the BE describes this event 
by an electron following the dashed line. This effective 
trajectory (from A to B) is longer than the real one (from 
A to B) by A s = 2VR 2 - b 2 . 

One can include the finite size of impurities into the 
transport equation in a manner to parallel traps. We 
approximate the trajectory of the electron by the effective 
trajectory AB. Since, following the a real trajectory AB, 
the electron reaches a next collision sooner by a time 
A s /u, we introduce into transport equation (^|) a negative 
time delay At — —A s /u. Here u is an electron velocity. 

For the hard-sphere impurities, transport equation (j^) 
with the negative collision delay is only an approxima- 
tion. Let us check how this approximation works for the 
correlated density. The classical scattering rate on hard 
spheres reads 

P P k^ [ -^-S(\k\-\p\)c'u^smd, (7) 

where c' = -ZVj mp /f2 is an impurity concentration (number 
of impurities N{ mp per volume £1), u = k/m is an electron 
velocity, and $ is a scattering angle, pk = \k\ \p\ cosd. The 
inverse lifetime follows from (^) as 



Here, we have denoted flu 



dp 



Ppk = c'wkB 2 . 



(8) 



The collision delay A t — — A s /u in angular coordinates 
reads 



\ t = --y/B 2 -b 2 = - -flsin-. 
u u 2 



(9) 



The correlated density from (^|) results 
dp 



(2tt) 3 
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/47T o 

-c' — B 3 n h 



imp 



Tiff 



(10) 



ric'^-B 3 



N imp fR 3 the 



total volume of impurities. 

The physical content of the correlated density can be 
demonstrated on the equation of state. The number of 
electrons which hit the surface of the sample is given by 
the density of freely traveling electrons nf re e, therefore 
the pressure P is given by the equation of state 



P = n llcc k B T. 



fill 



^From (g|) and ( JXOj ) we find relation between the free den- 
sity and the total number of particles N 



N 



(12) 



The equation of state flll| ) thus takes the form of the van 
der Waals equation 



P(fi - r> imp ) = Nk B T. 



(13) 



Briefly, the negative collision delay simulates for the ex- 
cluded volume in the van der Waals equation of state. 

Note that the correlated density (|10|) is negative. The 
density nf ree in the free space between impurities is higher 
than the physical density n what reflects that electrons 
are expelled from the volume of impurities. It is im- 
portant to distinguish which density (nf ree or n) is rele- 
vant for individual physical quantities. For instance, the 
charge density is given by n, while pressure relates to 

'*'£ rcc ■ 



III. COLLISION DELAY TIME IN QUANTUM 
MECHANICS 

The classical statistics shows that the non-locality of 
scattering events is approximatively described by the col- 
lision delay. This concept is easily transferred to the 
quantum mechanics, where the collision delay has already 
been introduced by Wignerla He used the maximum of 
wave packet to identify motion of an electron. Now we 
apply Wigner's approach to a neutral impurity to esti- 
mate a magnitude of virial corrections. 

The scattering of electron by a single impurity is de- 
scribed by the Schrodinger equatior 



(u-Ho- V)(V>in + V>out) = 0, 



(14) 



where ^in(^) = expifca; is an incoming plane wave, with 
r = (x, y, z), ip ou t is the outgoing part, Hq is Hamiltonian 
of unperturbed crystal and V is the impurity potential. 
The incoming plane wave has to be an eigen state of 
the crystal, (u — Ho)^- m — 0, thus the energy equals the 
kinetic energy of the incoming plane wave, u> = e/-. Then 
( |l4| ) simplifies as 



A formal solution of equation (f5h reads^ 



(15) 
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where 



Vw = G R (e k )T R (e k )^ u 
G » = ' 



(16) 



(17) 



uj - Ho + iO ' 

is the retarded Green's function of the host crystal, and 
T R = V + VG R T R , (18) 

is the T-matrix. 

As a model potential of the neutraLimpurity we use 
the one proposed by Koster and SlaterHij 

V = 10>v<0|, (19) 

where |0) is a single orbital at the impurity site. We 
will use the convention that lowercase denotes local el- 
ements of operators (that are in uppercase) throughout 
the paper. For the Koster-Slater potential, the T-matrix 
is also restricted to the selected orbital, T R = 10)^(01, 
and reads 



t 1 



v + v(0\G R \0)t R 



1 
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(20) 



To obtain the collision delay, we place the impurity in 
the initial of coordinates and express the wave function 
in the time representation 

Mr, t) = e***-"** - -^-t R (e k )e lkM - iekt . (21) 
2?r r 



We have used an asymptotic Green's function for large 
r, see Ref. EL 



(r\G R (e k )\Q)=- 



2?r M 



ik\r\ 



(22) 



to evaluate the outgoing wave from (|l^). This approx- 
imation holds for energies e k in the parabolic region of 
the band structure, e k — k 2 /2m. The first term in ( pl| ) is 
the incoming wave ip- m and the second one is the outgoing 
part Vout- 

To see the time delay, we take a linear combination of 
wave functions ip so that the incoming part ip- m forms a 
wave packet of a narrow momentum width k — > 0, 



dpe 



e ikx-ie k t exp 



lp-ky . . . 
— -ri — ^ipx-ztpt 



-(x-utf 



(23) 



where u = k/m is an electron velocity. This wave packet 
passes the initial of coordinates at t = 0. A corresponding 
outgoing wave ipout reads 




h- I , / i 01 " 




The outgoing wave passes the initial of coordinates with 
the collision delay 



A, = Im 



1 dt R 
t R ~d~J 



(25) 



The collision delay (]25|) depends only on the energy of 
electron. This is because the Koster-Slater impurity has 
a single scattering channel of the s-symmetry. For a gen- 
eral potential V, collision delay A t depends also on the 
scattering angle as the classical collision delay (^). 

The collision delay (|2^) is a quantum counterpart 
of the classical collision delay (^|). Following analogy 
between the quantum and classical approaches to the 
Boltzmann-like transport equations, we introduce the 
collision delay (25) into the scattering integral in exactly 
the same way as in the classical case. In other words, we 
expect the transport equation to be of form (Q), however, 
scattering rates P kp and collision delay A t are extracted 
from quantum collisions. 

The rate of scattering by Koster-Slater impurities of 
concentration c (probability that impurity occupy a site) 
follows from the Fermi golden rule as 



P P k =c\t R (e k )\ 2 2ir6(e k ~e p ). 



(26) 



This scattering rate does not depend on the scattering 
angle, thus it can be also expressed in terms of the life- 
time T 



p P k = \ 2 -^m-\k\), 



(27) 



where r is conveniently evaluated from the T-matrix 
1 



c(-2)Imt R (e fe ). 



(28) 



A. Estimate of virial corrections 

^From a scattering by a single impurity one can esti- 
mate magnitude of virial corrections. Using formula (|^) 
with the quantum scattering rate (|27]) and collision delay 
(p5|), one finds correlated density 



(27T) 



(29) 



The magnitude of virial corrections is thus measured by 
a ratio A t /r. 

Note that the collision delay is independent from the 
impurity concentration, while the lifetime is inversely 
proportional to the concentration. Accordingly, A t /r ~ 
c, i.e., magnitude of virial corrections is controlled by the 
impurity concentration c. To be specific, we will assume 
impurity concentrations ~ 10 -6 per site. 

Now we estimate A t /r for a model local Green's 
functionB 



4 



(0|G »|0) 
2 ( h b 3 



+ hz' 



1 -z 2 ) — (l + b 1 z + b 3 z 3 )^/T 



W 



(30) 



Here, W = 6 eV is a half-width of a conductivity band, 
and parameters b\ = 1.2 and 63 = — .4 serve to model 
the local density of state to a shape resembling III-V 
semiconductors, see Fig. ||. 

The collision delay is very sensitive to a value of the 
impurity potential v. Using (pp|), one can rearrange the 
collision delay (|25|) as 



A* = -Im 



dojt R 



= Im 



l-v(0\G§\0)' 



(31) 



Apparently, the collision delay will be long for potentials 
for which the denominator 1 — w(0|Gg'|0) goes to zero. 
For these values of potential v, the impurity behaves like 
a resonant level close to the conductivity band edge. 

For model function (|30|), the real part of the local 
Green's function at the band edge oj = equals to 
—0.185 1/eV. For potentials v < —5.4 eV, an impurity 
has a bound state. For v > —5.4 eV, there is a resonant 
level. In our calculations we use value v = —5.35 eV. 
In Fig g, the energy dependence of At evaluated from 
( |3i"| ) is compared with the lifetime r from ( f28[) . In Fig |] 
one can see that A t /r ~ 0.1, therefore appreciable virial 
corrections appear already for assumed concentration of 
resonant level c = 10~ 6 per site. 

The strong dependence of the collisional delay on the 
position of the resonant level leads to a strong depen- 
dence of virial corrections on the impurity potential, see 
Fig. ||. Such changes of the impurity potential can be 
achieved for instance by a hydrostatic pressure.u The im- 
purity concentration and the hydrostatic pressure can be 
thus used to control magnitude of virial corrections. 



IV. QUASIPARTICLE PICTURE 

^From analysis of the scattering by the Koster-Slater 
impurity, we have found that the largest virial corrections 
appear for resonant levels. Resonant levels, however, also 
result in large values of the T-matrix, as one can see in 

Fig. |. At the band edge Ret R 400 eV and lmt R ~ 

—30 eV. In particular, the real part of the T-matrix is 
large compared to potential v = —5.35 eV. For such large 
values of the T-matrix, the impurity scattering affects the 
electronic band structure. To take this effect into account 
we have to treat electrons as quasiparticles. 



A. Averaged T-matrix approximation 

The multiple scattering by impurities has been de- 
scribed in detail already in 1960's within Green's 



functions B In the averaged T-matrix approximation 
(ATA) that corresponds to our approximation of scat- 
tering rates, the self-energy equals the averaged value of 
the T-matrix, 

S ii = cT R (32) 

Since T R = \ r )'t R ( r )( r \t we can write the self-energy 
as 



(33) 



B. Energy renormalization 

The quasiparticle energy that describes propagation 
given by effective "Hamiltonian" Hq + cT reads 

Rea R (e k ). 



£k — (-k 



(34) 

The imaginary part of the self-energy provides the life 
time 



- = -2Jma R (e k ) 



(35) 



which is identical to the Fermi golden rule value (27). 

In Fig. [7] one can see that effect of impurities with the 
resonant level on the band structure is rather profound. 
The major effect is the overall shift of the band. This 
shift does not influence bulk properties of homogeneous 
crystals because it is compensated by a shift of chemical 
potential. 

The energy renormalization leads to quasiparticle cor- 
rections to velocity 



de k 
dk 



(36) 



Taking the momentum derivative from (|3J) one finds the 
renormalized velocity as 



k 

u = z — , 
m 



where z is the wave-function renormalization 
dRca R (uj) 



z(k) = 1 + 



(37) 



(38) 



With respect to transport properties, the velocity 
renormalization is the most important quasiparticle cor- 
rection as it determines drift of quasiparticles between 
collisions. The wave-function renormalization as a func- 
tion of the impurity potential is presented in Fig. |^. 
There is a striking similarity of the wave-function renor- 
malization and the magnitude of virial corrections. 

In the above discussion we have ignored selfconsistency. 
In Appendix A it is shown that for the weak scattering, 
- — > 0, the above formulas are identical to those result- 
ing from selfconsistent treatment. In Appendix B we also 
derive a formula that connects virial correction 1 + A t /r 
with quasiparticle renormalization z. This formula ex- 
plains similar magnitudes of these two corrections. 
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V. QUASIPARTICLE BOLTZMANN EQUATION 
WITH COLLISION DELAY 

Similarity of magnitudes of quasiparticle and virial cor- 
rections show that both corrections have to be included 
in the transport equation within the same accuracy. It 
is quite easy to guess such a transport equation. The 
quasiparticle renormalization affects the drift between 
collisions, therefore it enters the transport equation as 
a renormalization of velocity ( |37| ) . The virial corrections 
enter the scattering integrals like in (||). The transport 
equation that includes both corrections reads 



dt 



m dr 

7~P~ 



dr dk 
dp 



f 



5(\p\-\k\)f(p,r,t- A t ). (39) 



Although this equation has the classical form (g), its com- 
ponents z, t, and At are given by quantum-mechanical 
microscopic dynamics. One can also view ( |39| ) as a phe- 
nomenologic equation with momentum-dependent pa- 
rameters z, At and r. 

Beside the transport equation, one also needs relation 
of observables to distribution function /. From the equa- 
tion of continuity one finds that the physical density in- 
cludes only virial corrections, 



dk 



1 



(40) 



while the density of particle current has only quasiparti- 
cle corrections 



J = 



dk k 
z—f 



(27T) 

dk de 
(2tt) 3 dk 



(41) 



^From the energy conservation one finds that the energy 
density includes both corrections 



E = 



dk 



! + —)(£ + 0)/. 



(42) 



Both conservation laws are in A 
^From the set of equations 



endix C. 

42) one can evaluate 
properties of electron gas or liquid in a similar manner 
as one uses the BE to this end. To demonstrate such an 
application, in the next section we evaluate the dielectric 
function. 



VI. DIELECTRIC FUNCTION 

The virial corrections influence a response of the sys- 
tem to perturbations. The time non-locality of the scat- 
tering integral emerges in non-stationary processes. The 



simplest but important process is linear screening of ex- 
ternal field described by dielectric function K r . 

The virial corrections enter the dielectric function in 
two ways, from the transport equation ( |39| ) and from 
functional (Hoj) . To demonstrate both mechanisms, we 
evaluate k t from its definition. 

An electrostatic external potential 

Mr,t)=4>oe iqx - iu ' t (43) 
creates a perturbation in the electron density 

h(r,t) =fie lqr ~ lujt . (44) 
The perturbation in density creates Coulomb potential 



(r,t) = 4>e iqr ~ iuJt = JL-nji*-™* 
nq 



(45) 



that adds to the external one so that the internal field 
reads 



(46) 



90 

; 



Here, k is permitivity of the host crystal. From definition 

(47) 

one finds the dielectric function to be 

e n . . . 

K r = l-— -. (48) 



A. Perturbation of quasiparticle distribution 

To evaluate perturbation h of the physical density n, 
we have to find the linear perturbation of the quasipar- 
ticle distribution, 



f{k,r,t) = f{k)^- iut , 



(49) 



caused by potential d>. To this end we use linearized 
transport equation ( p9[ ) 



-ilu + iz h - ) J(k) - iqcp — — — 

m t I ok 



12ir 2 



dp 



t k 2 J (2tt) 3 



6(\ P \ ~ \k\)e^f( P ). 



(50) 



The momentum derivative of the equilibrium distribution 

fo(k) = fFD(ek) reads 



dfo(k) _ k df FD {e k ) 



<:)k 



m dek 



(51) 



The perturbation / depends only on the absolute value 
of momentum |fe| and the angle between momentum k 
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and wave vector q. We denote s — J%L-t and s' — T^rk, 
and integrate over the energy conserving 5 function so 
that the transport equation simplifies as 



m t J m o£\k\ 



l-e^ t ds'f(\k\,s'). (52) 



With abbreviations = u, \q\ = q, and skipping argu- 
ment |fc| in distributions, equation ( p0[ ) reads 



■ , ■ i ^ tt \ ■ A d f FD 
—iu) + isqu H — ] f[s) — isqucp — - — 

T 0£ 



1 



-e- A * j_ds'f(s'). (53) 



The angular dependence of the distribution is easily 
found from (|5£ 



/(*) 



where 



isqucj)^- + ± e iujAt F 
—iuj + isqu + - 



f = 2 y_/ s/(s) ' 



(54) 



(55) 



is an angle- averaged distribution. Integrating over s, one 
finds from ( |54] ) a condition for the angle-averaged distri- 
bution 

F=(l + (l-«JT)J)(f>?^-e i " A *JF, (56) 

ae 



where 



J = -i-ln ' " 



qu 



2qur \uj + - + qu 



(57) 



Since the BE holds only for slowly varying fields, we 
can linearize in A t , e lwAf ~ l+iujA t . The angle-averaged 
distribution from (B6|) then results 



F 



de 1 + (l + iojA t )J' 



(58) 



Now the perturbation of the quasiparticle distribution is 
fully determined by d5j) and (Est) . 



B. Perturbation of density 

Perturbation of the electron density is found from ( f4(i| ) 
dk 



h = 2 



■f(k) 



1 fOO / A 

dkk 2 F(k) 1 + — 



(59) 



The factor of two stands for sum over spins. 

For simplicity we assume the limit of low temperature 

^R^-S(e-E F ) = -^-d(k-k F ), (60) 
oe zk 

where one can easily integrate out the momentum 



mk 



1+ A^ l + (l-iuJT)J 



T J 1 + (l+iujA t )J 



(61) 



k—kf 



Using ( |6l| ) in ( [i"8| ) one directly obtains the dielectric func- 
tion. 



C. Long wave length limit 

Now we focus on long wave length limit, q — > 0. To 
evaluate this limit from (W\h we first rearrange (pTl) as 



J = -- 



1 



1 — iuiT 



1 



1 — IUIT 



1 - 



qur 



(62) 



3 \1 — iuiT / 

In the long wave length limit, the dielectric function reads 

1 I K.7T Z V T / 



1 — tLUT V ' Tt \ T t 



(63) 



where z, r and A t are values at the Fermi level. 

In the static case u> = 0, the dielectric function is of 

form K r = 1 + |f . From (|63| ) one finds that the Thomas- 
Fermi screening length l/q s is 



e 2 mkp 1 



A, 



(64) 



Part e gives standard Thomas-Fermi screening, fac- 
tors - and 1 + provide quasiparticle and virial correc- 
tions, respectively. As one can see from Fig. ^, quasiparti- 
cle and virial corrections nearly equal, therefore they mu- 
tually compensate in the Thomas-Fermi screening length 



e 2 mkF 



KIT 



(65) 



For homogeneous perturbations, q — 0, the dielectric 



function is of form n r = 1 — 
plasma frequency u> p results 



From (H), the 



(66) 



There is only the quasiparticle correction z. 
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Note that virial corrections to the Thomas-Fermi 
screening q s and plasma frequency io p appear in rather 
paradoxical way. While the static screening has virial 
corrections [due to n[f], Eq. (p9)], the plasma frequency 
describing non-stationary behavior has none. This is be- 
cause in the homogeneous case, q = 0, the virial cor- 
rections from n[f] and the scattering integral mutually 
cancel due to particle conservation law. 



E. dc conductivity 

Compensation of virial and quasiparticle corrections 
also appears for dc conductivity <7dc- Although this com- 
pensation is a direct consequence of the dielectric func- 
tion, we discuss it in detail for its experimental impor- 
tance. 

The conductivity relates to dielectric function 



D. Virial correction to Fermi momentum 

The Thomas-Fermi screening length (^J) and the 
plasma frequency ( |66| ) are expressed in terms of the Fermi 
momentum. Additional virial corrections to those quan- 
tities appear if one rewrite them in terms of physical 
density n. 

In general, the Fermi momentum is a parameter of the 
quasiparticle distribution /, therefore it is always related 
to the free density. For the parabolic band, the Fermi 
momentum results from (||) as 



k F = \J 3-7P 



o"dc = lim ~iuJK(K r — 1). 



(73) 



This known re lati on can be recovered from the equation 
of continuity (C2) that yields icon — iqj = 0, where j is 
flow of particles. The electric field F results from the 
electrostatic potential as eF = iqifi. The conductivity 
then reads 



ej e u> n 

Ode = = i — t- 

b q z q> 



(74) 



Comparing ( |74| ) with ( J48| ) one recovers (73). 

Sending q — > and uj — > one finds standard relax- 



^frcc 



(67) ation time formula with the quasiparticle correction 



For sufficiently low density n, the ratio A t /r changes a 
little from zero to the Fermi momentum. In our case, 
At/r changes by 7% for density n = 10 16 cm~ 3 . This 
weak dependence allows us to take approximation 



A«(fc) = A t (k F ) 
r{k) r{kF) 



(68) 



^From (^0|) one finds that the free and physical densities 
relate as 



Ufr, 



and the Fermi momentum reads 



k F 



3n 2 n 



1 



A: 



(69) 



(70) 



In terms of the physical density, the Thomas-Fermi 
screening length (p5J) regains corrections 



e 2 m / 37T 2 ?i 



(71) 



In opposite, the plasma frequency in terms of physical 
density regains its free-particle value 



2 

en 



Km 1 At 



(72) 



Cdc 



37r 2 m 



(75) 



In terms of the physical density, 



Ode 



e 2 nT 



m 1 



(76) 



virial and quasiparticle corrections mutually compensate. 



VII. SUMMARY 

We have shown that for scattering by resonant levels of 
neutral impurities the virial and quasiparticle corrections 
are of the same magnitude. We have proposed an intu- 
itive modification of the BE that includes both correc- 
tions. Proposed modification of the BE has quasiparticle 
corrections in the drift term (as in the Landau theory) 
and virial corrections in the scattering integral (as in the 
classical theory of dense gases) . The modified BE can be 
solved as simply as the standard BE. 

An interplay of virial and quasiparticle corrections has 
been discussed on the dielectric function. Various com- 
pensations of virial and quasiparticle corrections has been 
demonstrated on the static screening, plasma frequency 
and dc conductivity. Careful measurements of the dielec- 
tric function in III-V semiconductor with resonant levels 
can reveal this interplay. Sensitivity of resonant levels 
to a hydrostatic pressure makes possible to control mag- 
nitude of virial and quasiparticle corrections in a single 
sample. 
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The way we have introduced the transport equation 
([39|) does not guarantee its validity. There are two fun- 
damental questions one has to ask: (i) Has the time non- 
locality of the quantum-mechanical scattering integral re- 
ally the same form as the classical collision delay? (ii) Are 
quasiparticle and virial corrections included in consistent 
manner? No intuitive argument can give satisfactory an- 
swers to these questions. To prove yes-answers for both 
questions one has to recover the transport equation from 
quantum statistics. Such a microscopic theory is in the 
second paper of this sequence. 
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e k = £k + Recr R u (e k ). 



Using (A5), one can rewrite ( |A6| ) as 

e k = e k + Rea R (e k - Rea R u ) 



(A6) 



(A7) 



In argument of the self-energy we use definition (A6) to 
recover 



e k = e k + Rea R {e k ). 
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2. Wave- function renormalization 



Within the self-consistent treatment, the wave- 
function renormalization results asEj 



APPENDIX A: SELF-CONSISTENT ATA 

Here we show that for weak scattering, ^ — > 0, non- 
self-consistent formulas from Sec. VB result also from the 
self-consistent generalization of the self-energy. 

The effect of impurities on the band structure affects 
propagation of the electron in the region between im- 
purities. This change should be included in the propa- 
gator that enters the T-matrix. Instead of the unper- 
turbed propagator G R , the T-matrix in the self-energy 
(B2T) should be constructed from the full propagator G R 



Cit = v + V(0\G R \0)t R l{ , 



where G R is given by the Dyson equation 
r^R _ r*R _i_ nRs^R r*R 

Lr — U Q -+- (_r 2j self Lr , 



and 



'self 



ct 



self- 



(Al) 



(A2) 



(A3) 



This approximation is called the self-consistent ATA. El 

For the KS impurity, where the self-energy has no 
momentum dependence, the self-consistency is simply 
achieved by a complex shift of the energy argument 



G R (u J )=G R (u J -a R { ). 



(A4) 



,/From (|A4j) one finds that t R u (uj) =t R (uj — a R u ), there- 
fore 



a R { (u) = a R (u-a R l{ ). 



(A5) 



<9Rccr 1P 

1 self 



(A8) 



^From ( |A5| ) follows 



dRea R lf 



dRea 1 



dRea R lf 



which can be rewritten as 



1 



dRea 1 



dRca r 



(A9) 



(A10) 



The wave- function renormalization (AS) is thus identical 
to (H). 



3. Lifetime 



Within the salf-consistent treatment, the inverse life- 
time results aa!3 



- = zlma R u {e k ). 

T 

^From ( |A5| ) and ( |A6| ) one finds 
= <r R (e k - ilmcr^if) 



(All) 



a R (e k ) - Hm(x R u (e k ) 



da R 



duo 



(A12) 



In the last line we have used linear approximation in the 
imaginary part of the argument which holds for weak 
scattering, 7^0. 
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Imaginary part of (A12) 



ImoSfCefc) = Im<7 R (e fc ) - Ima R lf (e k ) 
can be rearranged as 



<9Rccr fl 



Formula ( All ) is thus identical to ( pq 



uj=e k 

(A13) 



(A14) 



for particle-particle interaction), it is quite natural that 
they do not include virial corrections. 

The scattering by resonant levels is far from the Born 
approximation, for model and parameters we cons ide r 
here \t | ~ 100 x \v\. In this case, the second term of (B3) 
is of the order of 10~ 3 while the first one is of the order 
of 10 _1 . Accordingly, the second term can be neglected, 
i.e., virial and quasiparticle corrections are of the same 
magnitude. 



APPENDIX C: CONSERVATION LAWS 



APPENDIX B: QUASIPARTICLE VERSUS 
VIRIAL CORRECTIONS 

Although virial and quasiparticle corrections describe 
different features of the quasiparticle transport, both of 
them are linked to energy derivatives of the T-matrix. 
From this link follows similarity of their magnitudes. 

One can rearrange formula (|2q ) in the way that reveals 
the relation of virial correction 1 + — to wave-function 



Im 



1 or 



renormalization z. Writing Q2q ) as At 
h, ~ Wl^Jj an< ^ t ^ ie mverse lifetime ( p8| ) as ^ 



2clmi 



ic(t R — i ), [t is complex conjugate to t R ], 



one obtains that 



ic(t R - t A ) 



2i 



1 dt R 
t R ~d^ 



1 dt A 



- —It « A\ c(t A dt R t R dt A 

~ &j2 (< +t >~ 2 \t R dto ~ t A du 



dRea 1 
duj 



c (t A dt R 
2 (t R ~dJ 



t R dt A 



(Bl) 



1. Equation of continuity 

In inhomogeneous and non-stationary system, there 
are currents j due to which local density of electrons 
n changes. Here we prove that density (|40j) and current 
(|Il]) are consistent with the BE ( p^ ) obeying the equation 
of continuity. 

Under integration over momentum, term |£ turns into 
drl Q t ° , and the scattering integrals turn into — ^| 
Sec. II. Equation (|39|) then yields 



see 



dn 
~dt 



dk 

Wr- 



it df d^df\ 

m or or ok I 



The second term in the brackets vanishes what follows 
from its integration by parts, 



dn 



dk k 

z—f = 0. 



c) 



dr J (2ir) 3 m 



(C2) 



The second term is the divergency of the current j given 
by (^l|). Equation (C2) is thus the equation of continuity. 



where we have used that Recr^ 
one finds 



-(t R +t A ). From (E 



r)t R f) 



which substituted into (pll) provides 



A* _ dRea R 
t du 



c\t R fRe^-(0\G R \0). 



(B2) 



(B3) 



Formula (B3) makes connection between virial and 
quasiparticle corrections. One can see that at least two 
limiting regimes can be distinguished according to rela- 
tive values of the first and the second terms in (B3). 

For weak potentials when the self-energy can be 
treated in the Born approximation t w v, i.e. a R = 
ctj 2 (0|Gq : |0), the virial corrections vanish because the 
first and the second terms mutually cancel. In contrast, 
the quasiparticle corrections remain. Since most of qua- 
siclassical transport equations have been derived within 
the Born approximation (or single-loop approximation 



2. Conservation of energy 

Here we prove that for homogeneous system the energy 
of electrons ( fl2|) changes with the field <f> in a consistent 



way, i.e., 
equation 



dE _ „d± 

■ dt 



First we take time derivative of 



dE 
~dt 



dk 8<j> f A t 
{2nf df'\ t 

(2^ (£ + ^i 



1 + ^ 

T 



(C3) 



The second term of (C3) vanishes. To show this, we 
multiply the BE (^) with the quasiparticle energy e + 4> 
and integrate over momentum 



dk 



} 8_l 

1 dt 



dk , ,.df A t 

W (e + *I7- (C4) 



We have used that contributions of non-gradient terms 
of the scattering integrals mutually cancel because of the 
energy-conserving S function. 
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In the first term of (C3) we take out the field, the rest of 
the integral is just density ( ft"o| ) . We have thus proved that 
the total energy changes in thermodynamically consistent 
way 



dE 
~dt 



at' 



(C5) 
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FIG. 1. Scheme of classical collision. The circle represents 
limits of an impurity potential, the full line is an electron 
trajectory, the dashed lines are extrapolations of incoming 
and outgoing trajectories. 



FIG. 2. Collision with a hard sphere. The full line is a real 
electron trajectory, the dashed lines are effective trajectories 
used within the BE. 




FIG. 3. Local Green's function. The imaginary part of 
local Green function (dotted line) has a low density of state 
at the band edge uj ~ corresponding to low effective mass 
m = 0.2, and high shoulder at higher energies simulating for 
satellite minima. The real part (full line) has finite value 
~ —0.185 1/eV and is nearly flat at the vicinity of the band 
edge. 
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to (meV) to (meV) 



FIG. 4. Collision delay and lifetime as function of energy. piG. 6. T-matrix near band edge. The real part Ret R is in 

Except for a shift by an order of magnitude, the collision f u u ii ne> the imaginary part lmt R is in dashed line, 
delay (full line) has nearly the same energy dependence as 
the lifetime (dashed line) . 




FIG. 5. Virial and quasiparticle corrections for electron en- 
ergy 100 meV above band edge and concentration c = 10~ 6 . 
Virial correction (dashed line) is greater than 1 for resonant 
levels, v > — 5.4 eV, what corresponds to positive collision 
delay. Quasiparticle renormalization z (dotted line) nearly 
equals the virial correction. In fact they differ less than by 
0.8%. 



FIG. 7. Quasiparticle energy (full line) for resonant levels, 
v = —5.35 eV, of concentration c = 10~ 6 . Bare kinetic energy 
tk (dotted line) serves as an eye guide. 
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